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Abstract — As the lung cancer is the leading cause of cancer 
death in the medical field, Computed Tomography (CT) scan 
of the thorax is widely applied in diagnoses for identifying 
the lung cancer. In this paper, a technique of rotation invariant 
with Local Binary Pattern (LBP) for segmentation of various 
lung nodules from the Lung CT cancer data sets is used. This 
is tested on various lung data sets from teaching files of 
Casimage database and National Cancer Institute (NCI) of 
National Biomedical Imaging Archive (NBIA). The results 
show the segmented nodules with clear boundaries, which is 
helpful in diagnosis of lung cancer. Further, the results are 
compared with the watershed segmentation method, which 
shows that LBP based method yields better segmentation 
accuracy. 

Index Terms — Image Segmentation, Lung Nodules, CT, 
LBP, MATITK. 

I. Introduction 

Cancer is the second most common cause of death in the 
United States after heart disease, accounting for 1 in every 4 
deaths [ 1 ] . As reported by American Cancer Society on can- 
cer, about 1 .6 million Americans were diagnosed with cancer 
in 2011 and 57 1,950 died of cancer, more than 1,500 peoplea 
day. In the past two decades numerous screening studies 
have been conducted worldwide to study early indications 
of lung cancer. Indeed, survival of lung cancer is strongly 
dependent on accurate and early diagnosis [2]. For identify- 
ing the lung nodules, Computed Tomography (CT) scan of 
the thorax is widely applied in diagnosis. The rapid develop- 
ments in chest CT acquisition techniques have been followed 
by a sharp increase in research on computer analysis of tho- 
racic CT scans. Compared to other modalities, CT excels in 
the imaging of the lungs [3] . The preprocessing step of most 
Computer-Aided Diagnosis (CAD) systems for identifying 
the lung diseases is lung segmentation [4] [5] . The purpose 
of the segmentation is to separate human body regions from 
background and make an initial classification showing the 
right and left lung clearly. For the detection of any lung dis- 
ease, it is further segmented to extract the exact ROI like lung 
nodules in the case of cancer detection [6] [7]. Hu et al. [8] 
described a method of global thresholding for the segmenta- 
tion. Pohle and Toennies [9] suggested adaptive region grow- 
ing method for segmentation of medical images. Most of seg- 
mentation methods are based on morphological operations 
[10]. But the results are not satisfied as judged by physi- 
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cians. Various automatic tools like Analyze 10.0 [12], Mazda 
[13], YaDiv [14] on CT are also suggested. Vincent Chu et al. 
[15] proposed MATLAB based tool called MATITK that uses 
the watershed algorithm for segmentation. This has given 
good results in lung segmentation especially for high vol- 
ume dataset [11]. There is no loss of lung nodules but prob- 
lem was the clear visualization. However, many papers are 
published on plain rotation invariance analysis [21,22,23]. A 
number of techniques incorporating invariance with respect 
to both spatial scale and rotation have been developed [24, 
25]. For the segmentation of Lung nodules on CT images, 
Rotation invariant Local Binary Pattern (LBP) [16] has shown 
very effective results. This helps in clear visualization of 
nodule boundaries which is important for doctors for analyz- 
ing the disease effects. 

II. Segmentation In Medical Imaging 

Segmentation has an important role in medical imaging as 
it helps in extracting the organ of interest. For the diagnosis 
of lungs, it is necessary to segment the chest images and 
extract the lungs in the preprocessing step and further the 
nodules are segmented with the different methods. One pul- 
monary CT image as shown in figure 1 is a precursor to most 
pulmonary image analysis applications. 




Figure l.The lung Cancer Image in which lung parenchymas seen as 
the dark region in the body. The surrounding tissue appears with a 
higher intensity 

A. Preprocessing Step: Lung Segmentation (to Separate 
Lung Parenchymas) 

By using the following steps, firstly lung area is separated 
from background [17] of the lung CT images as shown in 
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Figure 2. Separation of Lungs after the removal of mediastinum (a) 

Original DICOM CT Lung Image from LIDC (b) Thresholding 
point using histogram (c) Binary Image (d) Cleaning (e) Lung Mask 
(f) Lungs Extraction 

figure 2. 

1. Optimal Thresholding: The image is thresholded to 
separate low-density tissue (eg. Lung parenchymas) from 
fat. 

2. Background removal: The surrounding air, identified as 
low-density tissue, is removed. 

3. Cleaning: It is performed to fill the holes in the binary 
image. 

4. Lung mask: To extract lungs from background, lung mask 
is created. 

5. Lung Extraction: Output of step 4 is subtracted from the 
original image to provide separated lungs for further 
processing. 

B. Rotation Invariant Local Binary Pattern (LBP)for lung 
nodule Segmentation 

Local Binary Pattern (LBP) [16]: After finding out the lungs 
parenchymas from the preprocessing step, rotation invariant 
and gray scale invariance using Local Binary Operator is 
applied to get the lung nodules. The local binary pattern 
operator is theoretically simple, yet very powerful and gray- 
scale invariant method of analyzing textures [18] to do the 
segmentation. In practice, the LBP operator combines 
characteristics of statistical and structural texture analysis: it 
describes the texture with micro-primitives, often called 
textons, and their statistical placement rules [19]. Another 
important property is its computational simplicity, which 
makes it possible to analyze images in challenging real-time 
settings. 

Starting from the joint distribution of gray values of a 
circularly symmetric neighbor set of eight pixels in a 3x3 
neighborhood, there is an operator that is by definition 
invariant against any monotonic transformation of the gray 
scale. The dimensionality of the LBP feature distribution can 
be calculated according to the number of neighbors used. 
The definition of the LBP has been extended to arbitrary 
circular neighborhoods of the pixel to achieve multi-scale 
analysis and rotation invariance [16,20]. Rotation invariance 
is achieved by recognizing that gray scale invariant operator 
incorporates a fixed set of rotation invariant patterns. 
The value of the LBP code of a pixel (x c , y ) is given by: 
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LBP PJ} = r p -Jsig p _ gc WMx) ={J; %^ e: (1) 
Where 

g ; gray value of the center pixel 

g : gray values of the circularly symmetric neighborhood 

" g p (p=0,...,p-l). 
P: Image Pixels in the circle of radius R(R> 0) that form 

a circularly symmetric neighbor set 
2 P : binomial factor for each sign s(g - gj 

Since the LBP is invariant to monotonic changes in gray 
scale, it is supplemented by an orthogonal measure of local 
contrast. Below example shows how the contrast measure 
(C) is derived. The average of the gray levels below the center 
pixel is subtracted from that of the gray levels above (or 
equal to) the center pixel. Two-dimensional distributions of 
the LBP and local contrast measures are used as features. 
The operator is called LBP/C. 

An example of computing LBP and Contrast measure (C) 
in a 3*3 neighborhoods: 
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Pattern = 11110001 
LBP = 1 + 16 + 32 + 64 + 128 = 241 
C = (6+7+8+9+7)/5 -(5+2+l)/3 =4.7 
Important properties are: 

LBP is invariant to any monotonic grey level 

change 

Computational simplicity 
To remove the effect of rotation, i.e. to assign a unique 
identifier to each rotation invariant local binary pattern, 
rotation invariance is defined, formally: 

LBP PiE ri - mm {RaR(LBP P;R ,i)ji = P-l] (2) 

Where ROR(x, i) function performs a circular bit- wise right 
shift on the P-bit number x,i times. 

To quantify the varying performance of individual 
patterns attributes to the spatial structure of the patterns, a 
uniformity measure U ('pattern') is defined, which corresponds 
to the number of spatial transitions (bitwise 0/1 changes) in 
the 'pattern'. These patterns are called "uniform" because 
they have one thing in common: at most two one-to-zero or 
zero-to-one transitions in the circular binary code. However, 
only limited subsets of 'uniform' patterns are used instead of 
all rotation invariant patterns, which improve the rotation 
invariance considerably. So, this operator is called as LBP g 
R "" 2 . The use of 'uniform' patterns only is motivated by the 
reasoning that they tolerate rotation better because they 
contain fewer spatial transitions exposed to unwanted 
changes upon rotation. 

The LBP operator is an excellent measure of the spatial 
structure of local image texture, but it discards the other im- 
portant property of local image texture, contrast, since it de- 
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pends on the gray scale. If only rotation invariant texture 
analysis is desired, i.e., gray-scale invariance is not required. 
The performance of LBP is further enhanced by combining it 
with a rotation invariant variance measure VAR that charac- 
terizes the contrast of local image texture. 

Wm = ~£g£(M? - M) 2 .wktreii = (3) 

VAR g is invariant against shifts in gray scale. Since LBP and 
VAR are complementary, the joint distributions LBP s "" 2 /VAR g 
(aka LBPV) is very powerful rotation invariant measure of 
local image texture [16]. The combined operator LBPV 
separates the textures providing clear boundaries resulting 
in very efficient segmentation of image. 

III. Experimental Results 

In this study, the performance of rotation invariant LBP 
algorithm and watershed algorithm is analyzed on Lung Image 
Database Consortium (LIDC) public lung image database. 
The LIDC data have been collected from five different sites 
in the United States [28]. Watershed algorithm for volumetric 
segmentation using MATLAB Insight Segmentation and 
Registration Toolkit (MATITK) has used for lung 
segmentation [15, 17]. But still MATITK is under observation 
by some researchers as it is highly dependent upon the 
dataset. MATITK can be used only on 3D images, but Fig 3, 
shows the result of watershed segmentation on one slice of 
an image. There is no mechanism available to see all the 3d 
slices in MATITK. Using View3D' function, the segmentation 
of all the slices can be seen. 
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Figure 3. (a) Original Lung CT Image (b) Filtered Image (c) 
Segmented Image by Watershed Algorithm 

On the other hand, the rotation invariant LBP uniform pattern 
image LBP ' " and correspondingly histogram of the lung 
segmented nodules showing texture variance features after 
finding out the uniform pattern in figure 4(a) and 4(b) 
respectively. The final combined operator LBPV separates 
the textures with clear boundaries resulting in very efficient 
segmentation of image as shown in figure 5. 

Using SIENET Sky Dicom viewer tool [26] and RadiAnt 
DICOM viewer tool [27], segmentation of all the slices of 
medical images can be viewed efficiently In the Table 1 , the 
average accuracy 2 of the segmentation by the rotation 
invariant LBP Variance Operator is 98% over the accuracy of 
86% by MATITK method applied on 3 patient data sets of 
LIDC database [28]. 
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Figure 4. (a) Rotation Invariant LBP Uniform Pattern Image 
LBP"" 

(b) Histogram showing LBP Image using rotation invariant pattern 
of Image 




Figure 5. Segmentation of Lung Nodules using Combined LBPV 
operator 

Table I. Segmentation Average Accuracy 
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rV. Conclusion and Future Scope 

Though automatic segmentation is desirable but manual 
intervention in medical diagnosis is indispensable therefore, 
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Local Binary operator is used for gray-scale and rotation 
invariant texture classification. The approach is based on 
uniform patterns and nonparametric discrimination of sample 
and prototype distributions. Uniform patterns help to 
recognize texture microstructures such as edges. The operator 
performs well for the structural and stochastic patterns. Further, 
to show the efficacy of LBP based method, the results are 
compared with the watershed algorithm. The comparison 
indicates that LBP is very effective in segmenting the lung 
nodules on CT images and also gives the improved average 
segmentation accuracy over the watershed method which is 
further helpful in diagnosis of lung cancer. 

Segmented nodules can be further matched with differ- 
ent types of lung nodule templates for the retrieval of abnor- 
mal nodule, so that early retrieval of cancer can be done. 
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